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A  percolation  micro-model  is  developed  to  predict  the  effective  properties  of  the  composite  electrode 
consisting  of  poly-dispersed  electronic  and  ionic-conducting  particles.  The  analytical  expressions  for  the 
percolated  triple-phase-boundary  (TPB)  lengths,  hydraulic  pore  radius  and  intra/inter-particle  conduc¬ 
tivities  have  been  formed.  The  model  shows  that  the  percolated  TPB  lengths  of  a  composite  electrode 
consisting  of  poly-dispersed  particle  sizes  for  both  materials  (i.e.,  having  a  normal  distribution  with  a 
non-dimensional  standard  deviation  of  0.4)  are  about  32%  lower  than  that  of  the  composite  electrode 
with  mono-sized  particles.  And  the  higher  percolated  TPB  lengths  can  be  achieved  by  reducing  the  mean 
particle  size  of  each  material  and  narrowing  the  particle  size  distribution.  A  composite  cathode  may  bene¬ 
fit  from  it,  because  of  its  low  electrochemical  activity.  The  model  also  shows  that  the  composite  electrode 
with  large  mean  particle  size  radii  ratio  between  the  electrode-  and  electrolyte-materials  and  broad  par¬ 
ticle  size  distributions  can  not  only  provide  higher  inter-particle  ion  conductivity,  but  also  reduce  the 
percolation  volume  fraction  threshold  of  the  electrode-material.  And  this  may  be  helpful  for  a  composite 
anode  to  obtain  a  higher  overall  ion  conductivity. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  application  of  a  composite  electrode  with  both  electronic 
and  ionic  conducting  materials  has  been  regarded  as  an  attrac¬ 
tive  alterative  compared  to  the  electrode  with  single  material 
by  extending  the  electrochemical  reaction  zone  from  the  dense 
electrolyte  into  the  composite  electrode  [1-4].  This  in  turn  may 
significantly  increase  the  solid-oxide-fuel-cell  (SOFC)  performance 
[5-8].  A  typical  porous  composite  electrode  is  a  mixture  of  the 
electronic  conducting  materials  (i.e.,  lanthanum-doped  strontium 
manganate  (LSM)  for  cathode  and  Ni  for  anode)  and  ionic  con¬ 
ducting  materials  (i.e.,  yttria-stabilized  zirconia  (YSZ))  [3,5].  And 
a  proper  SOFC  composite  electrode  may  have,  but  is  not  limited 
to,  the  following  characteristics:  (i)  a  proper  thermal  conductiv¬ 
ity  that  matches  the  other  components,  (ii)  a  sufficient  mechanical 
strength  to  support  the  cell,  (iii)  a  percolated  porous  path  for 
the  reaction  species  to  transport  from  the  channel  into  the  reac¬ 
tion  sites.  The  ‘percolated’  in  here  is  defined  as  the  contiguous 
connection  through  the  whole  electrode  structure,  (iv)  a  per¬ 
colated  electronic  conducting  path  (through  Ni  or  LSM)  for  the 
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electrons  to  conduct  from  the  reaction  spots  to  the  current  col¬ 
lector,  (v)  a  percolated  ionic  conducting  path  (through  YSZ)  for 
the  ions  to  propagate  from  the  dense  electrolyte  to  the  reac¬ 
tion  spots,  (vi)  the  percolated  triple-phase-boundaries  (TPBs)  to 
support  the  electrochemical  reaction.  As  the  electrochemical  reac¬ 
tion  refers  to  electrons,  ions  and  reaction  species,  it  can  only 
occur  at  the  TPBs,  where  the  electronic  conducting,  ionic  con¬ 
ducting  and  porous  phases  coexist  (i.e.,  Ni-YSZ-pore)  and  each 
phase  belong  to  a  percolated  network,  respectively.  These  TPBs 
are  defined  as  the  percolated  TPBs  [9].  These  percolated  TPBs, 
however,  are  only  potential  electrochemical  active.  Flow  much 
part  of  them  can  be  active  again  rely  on  the  composite  electrode 
properties  and  the  thickness  of  the  composite  electrode  [4].  Con¬ 
clusively,  the  SOFC  performance  is  significantly  affected  by  the 
composite  electrode  properties.  This  in  turn  relies  on  the  partic¬ 
ular  microstructure  of  the  composite  electrode.  As  the  fabrication 
technology  of  a  composite  electrode  becomes  sophisticated,  the 
understanding  of  the  detailed  knowledge  of  microstructure  and 
its  effect  on  the  electrode  properties  are  receiving  considerable 
attentions  [5,10-13]. 

In  the  past  decades,  parametric  studies  for  the  composite  elec¬ 
trodes  have  been  widely  made  by  many  researchers  through 
experiments  [10,11,14-16],  random  packing  reconstruction  meth¬ 
ods  [17-20]  and  percolation  micro-models  [7,9,21-24].  (i)  The 
experimental  methods  rely  primarily  on  the  stereological  meth- 
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ods  [10,25-27].  Recently,  the  focused  ion-beam  scanning  electron 
microscopy  has  been  attempted  to  obtain  the  high  resolution 
3D  image  of  a  composite  electrode  [10,11,14,28-29].  All  of  them 
greatly  enhance  the  understanding  of  the  composite  electrode 
microstructure  and  are  helpful  for  developing  the  theory  model 
in  relating  the  microstructure  and  the  electrode  properties.  Cer¬ 
tainly,  the  theory  approaches  are  also  necessary,  because  the 
experimental  methods  are  still  expensive  and  time  consuming,  and 
require  the  use  of  the  hard-to-access  facilities  (i.e.,  focused  ion- 
beam,  synchrotron  radiation  facilities  and  laser  scanning  confocal 
microscopy)  [10,11,14,30].  (ii)The  random  packing  reconstruction 
methods  are  mainly  based  on  the  random  packing  of  sphere  par¬ 
ticles,  in  which  a  composite  electrode  is  considered  as  a  porous 
structure  formed  of  randomly  distributed  sphere  particles.  Ali 
et  al.  [19],  Schneider  et  al.  [18],  and  Golbert  et  al.  [17]  have 
reconstructed  the  binary  composite  structures  with  mono  size  par¬ 
ticles  for  each  material  to  evaluate  the  percolated  TPB  lengths. 
Recently,  Kenney  et  al.  [20]  have  further  extended  the  model  and 
reconstructed  a  binary  composite  structure  with  poly-disperse 
particles  to  predict  the  effective  electrode  properties.  However, 
only  limited  properties  have  been  calculated,  the  prediction  of 
some  other  electrode  properties  such  as  the  hydraulic  pore  radius, 
intra/inter-particle  conductivities  and  so  on  have  not  yet  been 
achieved  by  the  random  packing  reconstruction  methods  as  far 
as  we  know.  And  it  is  difficult  to  incorporate  the  reconstructed 
structures  into  the  cell-level  models  for  simultaneously  explor¬ 
ing  the  effects  of  various  microstructure  parameters  on  the  SOFC 
performance,  (iii)  The  percolation  micro-models  use  the  perco¬ 
lation  theory  and  coordination  number  theory  to  represent  the 
microstructure  of  a  composite  electrode  developed  by  the  random 
packing  reconstruction  methods.  And  analytical  expressions  are 
derived  to  relate  the  electrode  properties  and  the  microstructure 
parameters  [22,31].  Therefore,  the  percolation  micro-model  can 
be  considered  as  an  extension  of  the  random  packing  reconstruc¬ 
tion  method,  and  is  sufficiently  accurate  enough  in  predicting  the 
electrode  properties  from  the  microstructure  parameters  [22,24]. 
It  is  also  convenient  in  coupling  with  the  cell-level  models  [4]. 
Bouvard  et  al.  [32]  and  Suzuki  et  al.  [31]  have  developed  the 
micro-models  to  represent  the  microstructure  of  a  binary  ran¬ 
dom  packing  sphere  structure.  And  the  percolation  micro-models 
have  been  further  developed  by  Costamagna  et  al.  [22],  Chan 
et  al.  [21]  and  Janardhanan  et  al.  [24]  to  relate  the  composite 
electrode  microstructure  and  the  effective  properties.  Recently, 
Chen  et  al.  [9]  have  generalized  the  percolation  micro-model  of  a 
binary  random  packed  system  to  a  model  of  3-component  random 
packed  system  for  predicting  the  effective  properties  of  a  com¬ 
posite  electrode  with  LSM  particles,  coarse  and  fine  YSZ  particles. 
More  recently,  the  percolation  micro-models  have  been  further 
incorporated  into  the  cell-level  models  to  evaluate  the  effects  of 
electrode  microstructure  parameters  on  the  SOFC  performance 
[7,8,23]. 

In  the  early  work,  each  material  phase  of  the  composite  elec¬ 
trode  was  limited  to  the  mono-sized  particles.  But,  the  SOFC 
composite  ceramic  electrodes  are  typically  poly-disperse  powders 
that  each  material  may  have  a  broad  particle  size  distribution,  and 
this  can  have  significant  impacts  on  the  effective  electrode  prop¬ 
erties.  In  this  paper,  a  percolation  micro-model  is  developed  to 
predict  the  effective  properties  of  the  composite  electrodes  with 
a  distributed  particle  sizes  that  match  the  actual  ceramic  pow¬ 
ders.  Analytical  expressions  are  derived  to  predict  the  percolated 
TPB  lengths,  intra/inter-particle  conductivities,  etc.  based  on  the 
microstructure  parameters.  The  percolation  micro-model  is  vali¬ 
dated  through  comparing  with  the  random  packing  reconstruction 
method.  Furthermore,  all  the  effective  properties  are  represented  in 
non-dimensional  form  for  providing  generality  in  practical  appli¬ 
cations.  And  the  percolation  micro-model  is  particularly  focused 


on  the  applications  in  optimizing  the  SOFC  cathodes  and  anodes, 
respectively. 

2.  Theory  and  method 

2.1.  Percolation  micro-model  for  the  composite  electrodes  with 
multiple  particle  sizes 

2.1.1.  Coordination  number  and  microstructure  parameters 
Consider  a  multi-component  mixture  of  randomly  packed 

spherical  particles  with  M  type  particle  sizes.  The  average  number 
of  contacts  Zk  between  k-type  particles  and  all  neighboring  particles 
is  [9], 


M 


4 = y 

i=\ 


CD 


where  Zki  is  the  number  of  contacts  between  a  k-particle  and  all 
the  neighboring  ^-particles.  It  can  be  estimated  as  [9], 


zk,l  =  0.5 


,  ti/n 

rV  E  ti'/'k/rf 


(2) 


where  rk  is  the  radius  of  k-particle.  Z  is  the  average  coordina¬ 
tion  number  of  all  particles,  which  is  widely  assumed  to  be  6 
[21,22,33,34].  \j/k  is  the  volume  fraction  of  all  k-particles. 


2.1.2.  Effective  electrode  properties  and  coordination  number 

A  porous  composite  electrode  is  usually  assumed  to  have  two 
type  materials:  the  ionic-conducting  electrolyte-phase  material 
(labeled  as  ‘el’  such  as  YSZ)  and  electronic-conducting  electrode¬ 
phase  material  (labeled  as  ‘ed’  such  as  Ni  or  LSM).  Due  to  the 
ceramic  powders  are  typically  poly-disperse  powders,  a  composite 
electrode  consisting  of  m  type  particle  sizes  for  the  electronic- 
conducting  material  (designated  as  edi,  ed2, . . .,  edm)  and  n  type 
particle  sizes  for  the  ionic-conducting  material  (designated  as  eli, 
el2, . . .,  eln )  is  initially  considered  in  here  (m  +  n  =  M). 

Denoting  and  ^e\k  as  the  volume  fraction  of  edk-  and  e\k- 
particles,  respectively.  The  total  volume  fractions  of  the  electronic- 
conducting  and  ionic-conducting  materials  can  be  estimated  as, 

m  n 

V^ed  —  ^  ^V^edfc?  V^el  =  ^  (3) 

k= 1  k=  1 


And  the  relative  volume  fraction  of  each  type  particle  sizes 
within  the  electronic-conducting  and  ionic-conducting  materials 
can  be  represented  as 


M 


M 

fel  ' 


(4) 


The  percolated  TPB  lengths  per  unit  volume  (Z^pB  per)  is  pro¬ 
portional  to  the  contact  perimeter  between  ed/<-particle  and 
el^-particle  (/edfc  ei£ )» the  number  of  contact  points  per  unit  volume, 
which  is  the  product  of  the  numbers  of  ed^-particles  per  unit  vol¬ 
ume  (n^ )  and  their  coordination  numbers  with  other  el-particles 
(Ze dk,ei£),  and  the  probabilities  of  edk-  and  el^-particles  belonging 
to  the  percolated  clusters  (Ped/<  and  Pel£ ), 


\V 

ATPB,per 


el£nedkZedk,  e\£PedkPt 


el/* 


k= 1  1=1 


(5) 


As  indicated  in  Fig.  1 ,  the  contact  perimeter  (Zedfc  el  )  can  be  esti¬ 
mated  as  le dk,e\£  =  min(redfc,  re\  )  sin  0 ,  where  redfe  and  rel£  are 
the  radii  of  edk- particles  and  el^ -particles,  0  is  the  contact  angle. 
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Fig.  1.  Illustration  of  a  composite  electrode  with  both  coarse  and  fine  electrolyte- 
particle  sizes.  And  the  overlapping  between  edfc-particle  and  el^ -particle  is  described 
on  the  right  side. 


The  number  of  ed^-particles  per  unit  volume  within  the  entire 
composite  electrode  structure  can  be  estimated  as  [7,22], 


y  (1  0g)V^e  dk 

^  =  4^/3 


(6) 


where  <f>g  is  the  porosity. 

Since  only  the  percolated  TPBs  are  potential  electrochemi¬ 
cal  active,  it  is  the  percolation  probability  to  account  for  the 
effects  of  the  percolated  TPB  lengths.  For  a  binary  mixture 
with  two  mono  sized  particles,  edi-  and  eli -particles,  the  per¬ 
colation  probability  of  eli -particles  can  be  estimated  as  Pe\^  = 

[1  - ((3.764 -Zell,ell)/2)2-5]0'4  [22], 

However,  for  a  composite  electrode  with  n  type  particle  sizes  of 
electrolyte-material  as  indicated  in  Fig.  1,  the  percolation  proba¬ 
bility  of  elfc-particles  belonging  to  the  percolated  clusters  depends 
not  only  upon  the  coordination  numbers  among  the  same  size  of 
electrolyte-particles  (Zelfe  elfc),  but  also  on  the  mutual  coordination 
numbers  among  the  similar  materials  (Zelfc>  e\£  {k  i)).  In  order  to 
calculate  the  percolation  probability  of  el^-particles,  the  average 
number  of  contacts  between  el/<-particles  and  all  the  neighboring 
electrolyte-particles  (Zel/<el  =  Y^t= iZeik,ei£)  is  used  here, 


Peh  = 


1  - 


3-764  -ELi* 


'elfe.eli 


0.4 


(7) 


This  expression  is  validated  by  comparing  the  calculated  results 
through  the  percolation  micro-model  with  the  results  reported 
by  Kenney  et  al.  [20].  Similarly,  the  percolation  probability  of 
ed/<-particles  can  also  be  estimated  using  the  average  number  of 
contacts  between  ed^-particles  and  all  the  neighboring  electrode- 
particles  (Zed|(.ed). 

The  electrochemical  reaction  not  only  takes  place  within  the 
composite  electrode,  but  also  at  the  interface  of  the  composite  elec¬ 
trode  and  dense  electrolyte  [7].  And  the  percolated  TPB  lengths  per 
unit  electrolyte  surface  area  can  be  estimated  as, 

m 

ZTPB,per  =  ^2(27tredk  sin  9)nsedkPedk,  (8) 

k= 1 

where  the  number  of  ed/<-particles  per  unit  dense  electrolyte  sur- 
face  is  =  (1  -  <pg)i/edk2jrr^/3  [7,9], 

The  effective  intra-particle  ion  conductivity  based  on  the  geo¬ 
metric  dimensions  (<r^a,eff)  is  primarily  determined  by  the  effective 
relative  density  of  the  percolated  electrolyte-particles  (fei)  [35]. 
Jeon  et  al.  [7]  proposed  that  the  intra-particle  conductivity  can  be 


estimated  as, 


_tra,eff 

°el 

_tra,0 


_  tP  _ 

~~  5el  ~~ 


_k=  1 


(9) 


where  cr*;™’0  is  the  intrinsic  ion  conductivity.  The  Bruggeman  factor 
(/z)  is  used  to  include  the  effects  of  tortuous  conduction  phases. 
Similar  expression  can  also  be  applied  to  calculate  the  effective 
intra-particle  electron  conductivity  of  the  composite  electrode. 

For  a  electronic-conducting  electrode-phase  within  the  com¬ 
posite  electrode,  the  resistance  due  to  the  inter-particle  conduc¬ 
tivity  is  usually  negligibly  small  [9].  However,  the  inter-particle  ion 
conductivity  can  play  an  important  role  on  the  overall  ionic  conduc¬ 
tivity  for  a  ceramic  material  (i.e.,  YSZ),  especially  at  intermediate 
and  low  temperatures  [9,35].  And  the  particle  size,  inter-particle 
area,  doping  level  and  inter-particle  interface  thickness  can  signif¬ 
icantly  affect  the  inter-particle  ion  conductivity.  The  inter-particle 
ion  conductivity  (tx^r,eff)  may  be  evaluated  as  [9,21  ], 

Oreir’eff  _  Sk=lS€=l2aelk,el£(rel/<)2nelkZel/<,el£Pel^ 

1  j 

where  cr^r,°  is  the  intrinsic  ionic  conductivity  of  the  inter-particle 

interface.  aeifc, ei€  =  7r(min(reifc ,  re\e)  sin  0)2  is  the  surface  area  per 
contact  between  two  electrolyte-particles.  8  is  the  thickness  of  the 
inter-particle  interface. 

The  hydraulic  radius  is  used  mainly  for  describing  the  flow  in 
a  non-circular  tube,  and  is  also  an  essential  parameter  in  mod¬ 
eling  the  multi-component  porous  media  gas  transport  using  a 
Dusty-Gas  model  [36-38].  It  can  be  calculated  as  [9], 


where  x  is  an  adjustable  factor  that  represents  the  fraction  of  solid 
surface  areas  belonging  to  the  solid-pore  binary  phase  boundary 
areas  [9]. 


2.2.  Percolation  micro-model  for  the  uneven  distributed  particle 
sizes 


Actually,  for  a  practical  composite  electrode,  the  particle  sizes 
for  both  the  ionic-conducting  and  electronic-conducting  materi¬ 
als  may  have  uneven  probability  distributions.  In  present  paper, 
the  normal  distributions  are  used  to  represent  the  particle  size 
distributions  of  both  electronic-conducting  and  ionic-conducting 
materials.  Certainly,  many  different  kinds  of  probability  distribu¬ 
tions  can  also  be  applied  through  the  similar  approach. 

For  a  particular  material,  the  probability  density  function  of  a 
normal  distribution  with  a  specified  mean  radius  (r)  and  standard 
deviation  (cr)  is  given  as, 


f  = 


1 


exp 


-(r-r)2 
2  cr2 


(12) 


And  the  probability  of  particle  sizes  locating  in  the  range  of  r\ 
to  r2  can  be  calculated  as, 


-(r-r)2 

2o-2 


dr. 


(13) 


rr2  -|  rr  2 

p-l  mdr=^l, exp 

To  provide  considerable  generality  in  practical  application,  Eq. 
(13)  is  rewritten  as  a  function  of  the  non-dimensional  radius  (r/r) 
as, 

pr2/r  rr2l*  t 

p=f.  (14) 
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r/r 

Fig.  2.  The  probability  density  of  a  normal  distribution  as  a  function  of  r/r  with 
o'  =  0.2.  The  range  ( -  V2 o'  +  1 ,  V2 o'  +  1 )  is  discretized  into  9  type  particle  sizes  to 
represent  the  particle  size  distribution  of  each  material. 


13 


ved 


Fig.  3.  Comparison  between  the  percolation  micro-model  results  (lines)  and  the 
results  based  on  the  random  packing  reconstruction  method  (points)  [20].  A.^pB  per 
as  functions  of  \J/e d,  o'eA  and  o'e[  with  fed  =  rei  =  0.25  p.m,  6  =  29.5°  and  =  30  %. 


where  o’  =  o/f  is  defined  as  the  non-dimensional  standard  devia¬ 
tion.  In  this  way,  the  probability  density  of  a  normal  distribution  as 
a  function  of  r/r  with  &  =  0.2  is  shown  in  Fig.  2. 

As  illustrated  in  Fig.  2,  when  a  normal  distribution  is  used  to 
represent  the  particle  size  distribution  of  each  material,  there  is 
continuous  particle  size  distribution.  In  here,  only  the  range  from 
-\f2o'  +  1  to  \f2o'  +  1  is  considered,  because  of f'(-V2cr'  +  1)  = 
+  1)  =  (l/e)/'(l).  And  the  range  is  discretized  into  9  sec¬ 
tions,  therefore,  the  particle  size  distribution  of  each  material  is 
divided  into  9  type  particle  sizes  with  a  step  2^2o' /9. 

In  this  case,  the  probability  of  each  type  particle  size  ( pk )  can  be 
estimated  through  Eq.  ( 14).  Flowever,  it  is  necessary  to  mention  that 
the  cumulative  probability  density  of  all  9  type  particle  sizes  is  less 
than  1,  therefore,  the  probability  of  each  type  particle  sizes  should 
be  normalized  aspfc/^!=1p£.  In  Table  1,  the  non-dimensional  radii 
of  all  9  type  particle  sizes  and  their  normalized  probabilities  are 
estimated. 

And  the  relative  volume  fraction  of  each  type  particle  sizes  can 
be  estimated  as, 


<  = 


Pk[k 

£«=i  p*ri 


(15) 


Using  the  percolation  micro-model  described  above,  the  effec¬ 
tive  properties  of  a  composite  electrode  with  the  normal  particle 
size  distributions  can  be  predicted  conveniently.  There  are  only  two 
parameters  needed  to  specific  the  normal  distribution  function  of 
each  material’s  particle  sizes:  the  mean  particle  radius  (r)  and  the 
standard  deviation  (a)  (or  the  non-dimensional  standard  deviation 
(or'  =  a/r). 


3.  Results  and  discussion 

3.1.  Model  validation 

Fig.  3  illustrates  the  percolated  TPB  lengths  per  unit  volume 
(^tpb  per)  as  a  function  of  the  volume  fraction  of  the  electrode¬ 
material  (V^ed)  and  the  non-dimensional  standard  deviations  of 
the  electrode-  and  electrolyte-materials  (or'd  and  cr^)  with  red  = 
rel  =  0.25  pm,  0  =  29.5°  and  (pg  =  30%.  As  shown  in  this  figure, 
good  agreement  is  achieved  between  the  results  (lines)  based  on 
the  percolation  micro-model  and  the  results  (points)  based  on  the 
random  packing  reconstruction  method  reported  by  Kenney  et  al. 


[20].  It  demonstrates  the  suitability  of  the  proposed  percolation 
micro-model  in  relating  the  composite  electrode  properties  and 
the  microstructure  parameters.  It  should  be  pointed  out  that  there 
exist  several  transitions  in  the  curves  for  poly-disperse  case  (i.e., 
°ed  =  °ei  =  0,4  or  °ed  =  0*4,  °ei  =  reason  is  that  the  par¬ 

ticle  size  distribution  of  each  material  is  represented  by  9  type 
particle  sizes  and  each  type  of  particle  size  may  have  a  different 
percolation  volume  fraction  threshold  according  to  Eq.  (7).  There¬ 
fore,  each  type  of  particle  size  may  cause  a  small  transition  in  the 
curve.  These  transitions  can  be  eliminated  by  using  more  types  of 
particle  sizes  to  represent  the  particle  size  distribution. 

3.2.  Percolated  TPB  lengths 


To  provide  a  considerable  generality  of  the  percolation  micro¬ 
model  in  practical  application,  all  the  composite  electrode 
properties  are  presented  in  non-dimensional  form  basing  on  the 
underpinning  functions.  In  this  way,  the  results  are  not  specific  to 
a  particular  composite  electrode  structure,  and  once  the  physical 
parameters  are  given,  the  relevant  effective  properties  can  be  easily 
extracted. 

Substituting  the  contact  perimeter  and  particle  number  load¬ 
ing  into  Eq.  (5),  the  percolated  TPB  lengths  per  unit  volume  can  be 
rewritten  as, 


-EE! 


3  min(red. ,  rei  ) 


0  —  0g)  sin  0i/fedkZedk!  ei  PecjfcPei„ . 


(16) 


ed  k 


It  shows  that  ^pB  per  depends  upon  Zed(..  ei< .  Pedk  and  Pel( ,  which 
are  functions  of  the  volume  fraction  of  electrode-material  (V^d)* 
mean  particle  radii  ratio  of  the  electrode-  and  electrolyte-materials 
(fei/fed)  and  non-dimensional  standard  deviations  (c^d  and  <r^). 

The  non-dimensional  form  of  A^pB  per  can  be  written  as, 


TPB, per  ' 


ATPB,per 


(1  -<pg)  sin  6/r^ 


(17) 


Flere,  the  non-dimensional  percolated  TPB  lengths  per  unit  vol¬ 
ume  (X^pB  per)  is  independent  of  the  particular  particle  radii  (red  and 
rel),  contact  angle  (0)  and  porosity  (0g).  Flowever,  it  doses  depend 
upon  the  non-dimensional  standard  deviations  (cFd  and  cr^),  vol¬ 
ume  fraction  of  the  electrode-material  ( i^ed )  and  mean  particle  radii 
ratio  (re,/red). 
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Table  1 

The  non-dimensional  radii  of  9  type  particle  sizes  and  their  normalized  probabilities. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

rfr 

1  8^/2c^, 

1  9 

-i  6V2 a' 

1  9 

-1  4V2  o’ 

1  9 

-1  2V2 a’ 

1  9 

1 

-1  ,  2V2 & 

1  "T"  9 

-1  ,  4V2  a' 

1  "T"  9 

i  ,  Csfla' 

1  “T  9 

-1  1  8\/2 o' 

1  "T"  9 

P 

6.77% 

9.53% 

12.2% 

14.1% 

14.8% 

14.1% 

12.2% 

9.53% 

6.77% 

Fig.  4  shows  how  XjPB  per  depends  upon  ^ed,  cr'ed  and  or'  {  with 
rel/red  =  1-  As  indicated  in  Fig.  4,  the  maximum  X^pB  per  may  sig¬ 
nificantly  decrease  with  the  increasing  of  the  non-dimensional 
standard  deviation  of  each  material  (a'd  and  o'e{),  which  repre¬ 
sents  the  particle  size  distribution  of  each  material.  For  example, 
XjPB  per  with  the  particle  size  distributions  of  o'ed  =  o'el  =  0.4  and 

cr'ed  =  cr'j  =  0.6  are  32%  and  48%  lower  than  X^pB  per  with  mono 
particle  size  distributions  (or' d  =  or'  j  =  0)  with  the  same  mean  parti¬ 


cle  radii  and  x//ed  =  0.5,  respectively.  Similarly,  the  maximal  X^pB  per 
with  a'ed  =  0.4,  cr'j  =  0  is  21%  smaller  compared  to  the  maximal 
^TPB.per  with  =  cr^,  =  0. 

It  is  important  to  note  that  for  a  composite  electrode  with  red  = 
rel  and  cr'd  =  &v  the  volume  fraction  of  the  electrode-material 
(^ed)  for  the  maximum  X^pB  per  will  be  0.5.  For  the  case  red  =  rel 
but  cr'ed  =f=  cr'eV  however,  i/red  for  the  maximum  A^pB  per  may  shift. 


V^ed  for  the  maximum  XjPB  per  will  shift  from  0.5  to  a  higher  value 


for  cr'ed  >  a'eV  And  x//ed  for  the  maximum  X^pB  per  will  shift  from  0.5 
to  a  lower  value  for  o'ed  <  o'el. 

It  is  widely  agreed  that  there  is  a  threshold  in  the  vol¬ 
ume  fraction  of  the  electrode-  or  electrolyte-material,  percolated 
threshold  (i/Ad  or  t/Aj),  above  which  the  particles  form  a  continu¬ 
ous  network  through  the  whole  structure  [22].  Thus  the  volume 
fraction  of  the  electrode-material  should  be  chosen  within  the 
range  (i/Ad  <  V'ed  <  1  -  V^i)* In  ^is  case,  both  the  electrode-  and 
electrolyte-materials  form  the  connected  networks  through  the 
entire  composite  electrode  structure.  This  phenomenon  is  sup¬ 
ported  by  the  shapes  of  the  curves  shown  in  Fig.  4.  There  are 
volume  fraction  thresholds  for  both  the  electrode-  and  electrolyte- 
materials,  below  or  above  the  thresholds  (i/Ad  and  1  -  t/Aj)  there 
are  no  percolated  TPB  lengths.  And  it  can  also  be  found  that  the 
increasing  of  the  non-dimensional  standard  deviations  will  extend 
the  range  of  the  volume  fraction  thresholds  (i/Ad~(l  -  t/Aj)).  For 
example,  the  range  of  T/Ad~(l  -  V^i)  maY  extend  from  (0.3-0.7)  for 
(°ed  =  ^el  =  °)  t0  (°-2-°-8)  for  °ed  =  <1  =  0.4. 


Fig.  4.  The  effects  of  f  ed  and  non-dimensional  standard  deviations  (cr^  and  <r'j)  on 

^TPB,per  With  fei/fed  =  1. 


Fig.  5  illustrates  how  X^pB  per  depends  upon  \j/ed  and  rel/red 


with  o'ed  =  o'e j  =  0  and  o'ed  =  cr'j  =  0.4,  respectively.  As  indicated 
in  Fig.  5,  V'ed  for  the  maximum  X^pB  per  may  shift  from  a  higher 
value  to  a  lower  value  with  the  increasing  of  rel/red.  Fig.  5a  and 
b  also  shows  that  a  composite  electrode  with  the  broader  parti¬ 
cle  size  distributions  (represented  as  larger  o'ed  and  o')  may  have  a 
larger  percolation  volume  fraction  threshold  range  (i/Ad~(l  -  t/Aj)). 
For  example,  as  the  non-dimensional  standard  deviations  increase 
from  a’ed  =  <rd  =  0  to  a’eA  =  a’a  =  0.4,  the  range  of  i/^-P  -  fd)at 
rel/red  =  0-5,  1,  2,  4  may  vary  from  (0.46-0.82),  (0.30-0.70), 
(0.18-0.54),  (0.10-0.37)  to  (0.32-0.89),  (0.20-0.80),  (0.11-0.68), 
(0.06-0.51),  respectively.  And  the  maximum  X^pB  per  with  o'ed  = 
o'e j  =  0.4  are  25%,  32%,  25%  and  25%  lower  than  the  composite 


electrode  with  mono  particle  size  distributions  (cAd  =  cr'j  =  0)  at 
rel/red  =  0-5,  1,  2,  4,  respectively. 

Although  the  non-dimensional  \jPB  per  is  independent  of  the 
practical  particle  radius,  the  physical  A.^pB  per  actually  is  deter¬ 
mined  by  many  particular  microstructure  parameters.  For  example, 
let  0  =  29.5°  (j)g  =  0.3,  V^ed  =  0-5,  red  =  rel  =  0.25  p,m  and  o'ed  =  o'eX  = 
0.4,  Fig.  5b  can  be  used  to  calculate  A.^pB  per  =  11.6  x  1012  m-2. 
Therefore,  the  curves  for  non-dimensional  properties  are  very  gen¬ 
eral  for  practical  application,  and  the  relevant  effective  properties 
can  be  easily  extracted  as  physical  parameters  are  specified. 

Conclusively,  a  higher  X^pB  per  may  be  expected  for  the  compos¬ 
ite  electrode  with  small  mean  particle  sizes  of  both  electrode-  and 
electrolyte-materials  (red  and  rel)  and  narrow  particle  size  distribu¬ 


ted 


Fig.  5.  X^pB  per  as  functions  of  fed  and  fei/fed  with  (a)  =  0,  (b)  <x;d  =  a'el  = 

0.4. 
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Fig.  6.  The  non-dimensional  effective  intra-particle  conductivities  as  functions  of 
\fsed,  G'ed  and  &  with  rei/red  =  1.  The  upper  panel  for  d^a’eff  and  the  lower  panel  for 


tions  (<y  d  and  or'j).  And  this  conclusion  may  be  useful  for  the  SOFC 
composite  cathode  (i.eM  LSM  and  YSZ).  As  there  are  low  exchange 
transfer  current  density  based  on  per  unit  TPB  length  (low  electro¬ 
chemical  activity),  a  higher  X^pB  per  is  essential  for  enhancing  the 
cathode  electrochemical  performance.  Thus,  a  composite  cathode 
with  small  mean  particle  sizes  and  narrow  particle  size  distribu¬ 
tions  may  lead  to  the  optimized  structure. 

3.3.  Effective  intra-particle  conductivity 


The  intra-particle  conductivity  may  become  an  important 
parameter  for  the  charge-transfer  processes  within  the  composite 
electrode.  In  terms  of  the  TPB  lengths,  there  is  benefit  in  repre¬ 
senting  the  results  in  non-dimensional  form.  According  to  Eq.  (9), 
the  non-dimensional  intra-particle  electron  conductivity  may  be 
written  as, 


~tra,eff 
(7  A 
ed 


_tra,eff 

^ed 


^d’° 


(1  -4>gf 


(18) 


An  analogous  definition  can  also  be  used  for  the  intra-particle 
ion  conductivity  through  the  electrolyte-particles  (cf^a,eff). 

Fig.  6  illustrates  o^da,eff  and  df^a,eff  as  functions  of  the  vol¬ 
ume  fractions  of  the  electrode-material  (V^ed)  and  non-dimensional 
standard  deviations  [&ed  and  cr',)  with  rel/red  =  1.  As  indicated  in 

Fig.  6a,  cffda’eff  may  increase  with  the  increasing  of  \f/ed.  And  there  is 
also  percolation  volume  fraction  threshold  limit  of  the  electrode¬ 
material  (V^d)’  below  which  o^da,eff  is  zero. 

The  present  of  may  be  useful  for  the  composite  anodes 
with  Ni  and  YSZ.  As  there  are  high  electrochemical  activity  and 
high  specific  electron  conductivity  of  Ni  (o^da,°),  a  higher  effective 
ion  conductivity  is  essential  to  enhance  the  ions  transfer  capabil¬ 
ity.  This  in  turn  may  extend  the  electrochemical  active  zone  and 


decrease  the  anodic  potential  loss.  As  cf^a,eff  may  increase  with 
the  decreasing  of  \jred  (shown  in  Fig.  6b),  the  volume  fraction  of 
electrode-material  (^ed)  around  i/d,d  is  always  adopted  for  a  com¬ 
posite  anode  to  achieve  a  high  ionic  conductivity.  In  this  case,  both 
the  electrode-  and  electrolyte-particles  should  be  percolated  and  a 
high  effective  ionic  conductivity  can  be  achieved. 

The  effects  of  or'd  and  cr'e[  on  o^da,eff  and  a^a,eff  are  also  shown  in 
Fig.  6a  and  b,  respectively.  It  is  apparent  that  there  are  weak  effects 
of  standard  deviations  on  <5^da,eff  and  a^a,eff  at  a  high  ^red,  even 
with  the  broad  particle  size  distributions  ( or'd  =  =  0.6).  FIow- 

ever,  as  can  be  seen  from  Fig.  6a,  the  percolation  volume  fraction 
threshold  of  the  electrode-material  (^d)  ma Y  decrease  with  the 
increasing  of  the  non-dimensional  standard  deviations  (or'd  and 
&  x).  For  example,  when  the  non-dimensional  standard  deviations 
increase  from  cr'ed  =  cr'el  =  0  to  cr'd  =  cr'e{  =  0.6,  V^d  may  shift  from 
0.3  to  0.1 8.  In  this  case,  V^ed  around  ^d  =  0. 1 8  (for  cr'ed  =  cr'el  =  0.6) 
can  be  adopted  instead  of  V^d  =  0.3  (for  cx'ed  =  & x  —  0),  and  a  higher 
o^[a,eff  can  be  achieved.  This  in  turn  may  increase  cr^a,eff  from  0.59 
to  0.72,  approximately. 

Analogously,  there  is  also  volume  fraction  threshold  limit  in 
Fig.  6b  (1  -  t/^),  above  which  there  is  no  o^[a’eff.  And  1  -  ^  may 
shift  from  0.7  to  0.82  as  the  non-dimensional  standard  deviations 
increase  from  cr'd  =  a'eX  =  0  to  cr'd  =  a'e]  =  0.6. 

Fig.  7  shows  or*;da’eff  and  d^a,eff  as  functions  of  ^ed  and  rel/red 
with  <y  d  =  a'el  =  0.4.  As  illustrated  in  Fig.  7a,  for  a  given  \J/ed,  ^da,eff 
may  increase  with  the  rise  of  the  radial  ratio  rel/red.  And  o^a,eff, 
however,  may  decrease  with  the  increasing  of  rel/red  as  shown  in 
Fig.  7b. 

For  a  SOFC  composite  anode  with  Ni  and  YSZ,  a  volume  fraction 
of  the  electrode-material  ( ifed )  around  the  percolation  volume  frac¬ 
tion  threshold  \j/^d  is  always  chosen  to  achieve  a  high  intra-particle 


Fig.  7.  The  non-dimensional  effective  intra-particle  conductivities  as  functions  of 
yred  and  rel/red  with  cr^d  =  <f  =  0.4.  The  upper  panel  for  a^ja,eff  and  the  lower  panel 
for  cFja,eff. 
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Fig.  8.  The  non-dimensional  inter-particle  ion  conductivity  as  a  function  of  fed,  o'  Fig* 9*  The  non-dimensional  inter-particle  ion  conductivity  as  functions  of  f  ed  and 
and  ^  with  fel/fed  =  1.  Wred  with  ^  ^  =  0.4. 


ion  conductivity.  Although  cf^a,eff  may  decrease  with  the  rise  of 
fel/red  at  a  given  \j/ed  (shown  in  Fig.  7b),  i/Ad  may  also  decrease 
with  the  increasing  of  rel/red  (shown  in  Fig.  7a).  Therefore,  a  high 
^tra,eff  can  benefit  from  a  large  rel/red,  because  of  the  decreas¬ 
ing  of  i/Ad.  As  shown  in  Fig.  7a,  when  rel/red  increase  from  1  to  4,  a 
higher  o^a,eff  can  be  obtained  by  using  a  \j/ed  around  =  0.1  (for 
fei/fed  =  4)  instead  of  V^d  =  0.22  (for  rel/red  =  1).  Under  the  cir¬ 
cumstances,  cf^a,eff  may  increases  from  0.7  to  0.88,  approximately 
(shown  in  Fig.  7b). 

Based  upon  the  earlier  literature,  Chan  et  al.  [21  ]  reported  that 
the  specific  intra-particle  electron  and  ion  conductivities  for  LSM 
and  YSZ  are  around  cr*;da’0  —  10,000  S  m-1  (specific  resistivity 
around  0.0001  £2m)  and  cr^a,°  =  6.7  S  m-1  (specific  resistivity 
around  0.1 5  £2m)  at  945  °C,  respectively.  Considering  a  situation 
where  0g  =  O.3,  V^ed  =  0.5,  red/rel  =  1  \±m  and  cr'ed  =  cr'e{  =  0.4,  the 
effective  intra-particle  electron  and  ion  conductivities  can  be  eval¬ 
uated  as  cTgda,eff  =  1953  S  m_1  and  cr^a,eff  =  1.3  S  nrr1  through 
Fig.  7a  and  b,  respectively. 


3.4.  Effective  inter-particle  conductivity 


Fig.  8  shows  the  non-dimensional  inter-particle  ion  conductivity 
(d^r,eff)  as  a  function  of  \j/ed,  cr'ed  and  cr^  with  red/rel  =  1.  cf^r,eff  is 
defined  as, 


ter,eff 

^ter,eff  _  _ ^el _ 

el  0^’°^ dsin26»(l -(/>g)/5 


(19) 


For  YSZ,  the  inter-particle  thicknesses  8  have  been  estimated  to 
be  around  5  nm  [39]. 

As  illustrated  in  Fig.  8,  <5^r,eff  may  decrease  with  the  increasing 
of  the  volume  fraction  of  electrode-material  (V^ed)  with  fed/rel  = 
1.  For  the  narrow  particle  size  distributions  (i.e.,  a'ed  =  a'el  <  0.4), 
there  are  weak  effects  of  the  non-dimensional  standard  deviations 
on  <r^r,eff.  For  the  broad  particle  size  distributions  (i.e.,  &ed  =  cr'el  = 
0.6),  however,  the  effects  of  non-dimensional  standard  deviations 
on  o^r,eff  may  not  be  ignored  in  practice.  And  there  exist  about 
approximate  10%  increase  compared  to  the  composite  electrode 
with  mono  particle  size  distributions  (<y  d  =  <j'el  =  0). 

Fig.  9  shows  o^r,eff  as  functions  of  ^ed  and  fel/red  with  or'd  = 
j  =  0.4.  As  indicated  in  Fig.  9,  cr^r,eff  may  increase  with  the  rise  of 
fei/fed  for  a  given  ^ed.  Furthermore,  although  the  function  form 


of  the  non-dimensional  d^r,eff  is  independent  of  the  particular 
particle  radii,  the  physical  inter-particle  ion  conductivity  (T^r,eff  is 
proportional  to  the  mean  particle  size  of  each  material.  And  the 
physical  cr^r,eff  may  significantly  increase  with  the  rise  of  the  mean 
particle  size  of  the  electrolyte-material  (rel). 

Chan  et  al.  [21]  also  reported  that  the  specific  inter-particle 
ion  conductivity  for  YSZ  is  around  <x^r’°  =  0.05  S  m_1  at  945  °C 
(i.e.,  specific  resistivity  around  20  tl  m).  Consider  an  example 
of  0  =  29.5°  0g  =  O.3,  \//ed  =  0.5,  red  =  rel  =  0.25  pan  and  cr'ed  = 

=  0.4,  the  inter-particle  ion  conductivity  can  be  evaluated  as 

a“r  eff  =  0.87  S  m-1  by  Fig.  9.  Thus  <x^r’eff  =  0.87  S  nr1  has 
the  same  order  of  magnitude  with  the  intra-particle  ion  conductiv¬ 
ity  (<r^a,eff  =  1.3  S  nrr 1 ),  which  illustrates  that  the  inter-particle 
ion  conductivity  of  YSZ  can  play  an  important  role  on  overall  ionic 
conductivity. 

In  conclusion,  although  the  large  rel/red  and  non-dimensional 
standard  deviations  may  lead  to  a  low  X^pB  per,  they  are  helpful  for 
achieving  a  high  effective  ionic  conductivity.  A  composite  electrode 
with  a  large  rel/red  and  broad  particle  size  distributions  can  not 
only  provide  the  higher  (T^r,eff,  but  also  reduce  the  percolation  vol¬ 
ume  fraction  threshold  of  electrode-material  (V^d). As  a  comPosite 
anode  is  of  the  high  electrochemical  activity  and  intrinsic  electronic 
conductivity,  \/fe d  around  i/d,d  is  usually  used  to  achieve  the  higher 
intra-particle  ion  conductivity.  Therefore,  both  <7^a,eff  and  <7^r,eff 
may  benefit  from  large  rel/red,  a'd  and  a'eV  due  to  the  increase  of 
(T^r,eff  and  reduction  of  ^d-  Furthermore,  the  results  are  consis¬ 
tent  to  Itoh  et  al.  and  Kim  et  al.  finding  that  the  composite  anode 
using  larger  mean  particle  radius  of  electrolyte-material  and  broad 
particle  size  distribution  may  effectively  prevent  the  agglomera¬ 
tion  of  Ni  particles  and  increase  the  durability  of  the  SOFC  during 
long-term  running  [40,41  ]. 

4.  Summary  and  conclusions 

A  percolation  micro-model  has  been  developed  to  predict  the 
effective  properties  of  a  SOFC  composite  electrode  with  poly- 
disperse  particles.  The  properties  include  the  intra/inter-particle 
ion  conductivities  (<r^a,eff  and  cr^r,eff),  percolated  TPB  lengths  per 
unit  volume  (A^pB  ),  percolated  TPB  lengths  per  dense  electrolyte 
surface  (A^pB  per)  and  hydraulic  pore  radius  (rg).  These  properties 
depend  upon  the  physical  characteristics  of  the  microstructure, 
such  as  the  mean  particle  sizes  of  the  electrode-  and  electrolyte- 
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materials  (red  and  rel),  volume  fraction  of  the  electrode-material 
( V^ed)*  non-dimensional  standard  deviation  of  each  material  phase 
(<^d  and  cr'j),  intrinsic  material  properties  and  the  porosity  (0g). 
And  the  proposed  percolation  micro-model  can  be  conveniently 
incorporated  with  the  cell-level  model  for  optimization  of  the  com¬ 
posite  electrode  structure. 

The  percolation  micro-model  has  been  validated  through 
comparing  with  the  simulated  results  by  the  random  packing 
reconstruction  method  [20],  which  demonstrates  the  accuracy  of 
the  percolation  micro-model  in  dealing  with  the  composite  elec¬ 
trode  with  poly-disperse  particle  size  distributions. 

The  results  showed  that  a  higher  X^pB  per  may  be  expected  for 
the  composite  electrode  with  small  mean  particle  sizes  of  the 
electrode-  and  electrolyte-materials  (red  and  rel)  and  narrow  parti¬ 
cle  size  distributions.  And  a  composite  cathode  with  LSM  and  YSZ 
may  benefit  from  it  for  pursuing  a  high  SOFC  performance. 

The  results  also  indicated  that  a  composite  anode  with  large 
rel/red  and  broad  particle  size  distributions  can  not  only  pre¬ 
vent  the  agglomeration  of  the  Ni  particles,  but  also  increase  the 
inter-particle  ion  conductivity  and  reduce  the  percolation  volume 
fraction  threshold  of  the  electrode-material  (^d).  And  this  con¬ 
clusion  may  be  useful  for  a  composite  anode  with  Ni  and  YSZ  to 
achieve  a  high  overall  ionic  conductivity. 
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